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Background. An increasing proportion of people are living with multiple health conditions, or 
'multimorbidity'. Measures of multimorbidity are useful in studies of interventions in primary care 
to take account of confounding due to differences in case-mix. 

Objectives. Assess the predictive validity of commonly used measures of multimorbidity in rela- 
tion to a health outcome (mortality) and a measure of health service utilization (consultation rate). 

Methods. We included 95372 patients registered on 1 April 2005 at 174 English general practices 
included in the General Practice Research Database. Using regression models we compared the 
explanatory power of six measures of multimorbidity: count of chronic diseases from the Quality 
and Outcomes Framework (QOF); Charlson index; count of prescribed drugs; three measures from 
the John Hopkins ACG software [Expanded Diagnosis Clusters count (EDCs), Adjusted Clinical 
Groups (ACGs), Resource Utilisation Bands (RUBs)]. 

Results. A model containing demographics and GP practice alone explained 22% of the uncer- 
tainty in consultation rates. The number of prescribed drugs, ACG category, EDC count, RUB cat- 
egory, QOF disease count, or Charlson index increased this to 42%, 37%, 36%, 35%, 30%, and 26%, 
respectively. Measures of multimorbidity made little difference to the fit of a model predicting 
3-year mortality. Nonetheless, Charlson index score was the best performing measure, followed 
by the number of prescribed drugs. 

Conclusion. The number of prescribed drugs is the most powerful measure for predicting future 
consultations and the second most powerful measure for predicting mortality. It may have poten- 
tial as a simple proxy measure of multimorbidity in primary care. 

Keywords. Comorbidity, family practice, mortality, outcome assessment-health care, risk 
adjustment. 



Introduction 

The focus of primary health care in developed coun- 
tries now largely relates to supporting the manage- 
ment of long-term conditions. As the population ages, 
an increasing proportion of people are living with 
multiple health conditions, or 'multimorbidity'. There 
is a growing body of evidence about the epidemiology 
and consequences of multimorbidity, and a number 
of different measures of multimorbidity have been 
used in these studies. 1-4 Measures of multimorbidity 
are useful in both epidemiological and experimental 
studies of interventions in primary care, in order to 
take account of confounding due to differences in 
case-mix. 



Multimorbidity can be conceptualized and measured 
in several different ways. Some measures are based on 
simple counts of patients with multiple health conditions 
(which raises questions about the range of conditions or 
diseases that should be included), while other measures 
differentially weight diseases or use other approaches 
to take account of the burden of illness on an indi- 
vidual. 56 Most commonly used measures of multimor- 
bidity were originally developed and validated among 
selected patients in hospital. 6,7 The measures were also 
often developed to demonstrate relationships with a 
specific outcome, such as mortality, 8 but have since been 
used in studies with very different outcomes. 6 It can- 
not be assumed that measures of multimorbidity will 
be equally valid in populations encountered in primary 
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care settings, or are able to predict outcomes other than 
those for which they were developed. 

Several previous studies have compared the 
performance of different measures of multimorbidity 
to predict a specific outcome (such as mortality) in a 
specific population (such as people with hypertension). 
However, few previous studies have directly 
compared the performance of different measures of 
multimorbidity in relation to more than one outcome 
in a general primary care or community-dwelling 
population. These studies have been restricted to elderly 
populations 910 except for one study that explored the 
ability of five measures (the Charlson index, Elixhauser 
index, Chronic Disease Score, number of diagnoses 
and number of drugs prescribed) to predict mortality 
and hospitilization in adults living in Saskatchewan, 
Canada. 11 This study was based on linked administrative 
data such as billing and prescribing records. 

The purpose of this paper is to assess the predictive 
validity in primary care of commonly used measures 
of multimorbidity in relation to two different types of 
outcome: a health outcome (mortality) and a measure 
of health service utilization (primary care consultation 
rate). We also wanted to assess whether one measure 
provided greatest explanatory power in relation to both 
outcomes or whether different measures were needed 
for different outcomes. 

Methods 

Sample 

This study was based on routine anonymised electronic 
medical records from patients registered with prac- 
tices contributing to the General Practice Research 
Database (GPRD), which is now incorporated in the 
Clinical Practice Research Datalink (http://www.cprd. 
com). The GPRD contains the primary care records 
of around 5 million patients and is considered broadly 
representative of the general population in the UK. 12 
We used a stratified random sample of 99 997 individu- 
als aged 18 years and over and registered at one of 182 
practices on 1 April 2005 and obtained data on all con- 
sultations, diagnoses and prescriptions until 31 March 
2008. The sample was stratified by age, gender and 
practice. After excluding eight practices in which data 
about patient deprivation were unavailable, 95 372 
patients from 174 practices were used for analysis. 

Measures of multimorbidity 

Multimorbidity was operationalized in six ways. The 
first three multimorbidity indices discussed here used 
all historic diagnoses up until 1 April 2005. 

First, we constructed a simple count of 17 diseases 
included in the Quality and Outcomes Framework 
(QOF) 13 that have been classified as chronic in previ- 
ous research. 1 To identify whether an individual had a 



relevant disease we employed Version 11 of the QOF 
Business Rules. 14 

Second, we used the Charlson index score, in the recent 
adaption for use with Read coded data. 15 The Charlson 
index is a weighted score that includes 17 chronic diseases, 
where the weighting associated with each disease reflects 
the strength of its association with inpatient mortality. 8 

The next three multimorbidity measures were 
constructed using output from the John Hopkins 
University ACG Case-Mix System, which was devel- 
oped using administrative data in the USA. 16 The soft- 
ware allows the user to input Read Version 2 diagnostic 
codes. As our third measure of multimorbidity we used 
a count of Expanded Diagnosis Clusters (EDCs), which 
are groupings of diagnostic codes based on disease or 
organ area. The software identifies 264 separate EDCs, 
of which 114 were classified as chronic during work pre- 
viously described. 1 We counted the number of chronic 
EDCs ever recorded for each individual. 

Our fourth multimorbidity measure used Adjusted 
Clinical Groups (ACGs), which are mutually exclusive 
categories defined by age, gender and combination of 
Aggregated Diagnosis Groups (ADGs). ADGs are 
groupings of diagnostic codes based on clinical dimen- 
sions that determine the subsequent expected resource 
use. Examples of the dimensions used to classify diagno- 
ses into the ADGs are severity of the condition, expected 
duration of the condition, and the need for specialty care 
involvement. The ACG classifications used diagnoses 
recorded over a 1-year period (1 April 2004 to 31 March 
2005), as suggested by the software manual. 16 The age 
range of our sample meant that a total of 68 ACG cat- 
egories were populated (although the version of the soft- 
ware used identified up to 82 default categories, some of 
these were age specific). The fifth measure of multimor- 
bidity used Resource Utilisation Bands (RUBs), which 
are an aggregation of ACGs into six ordinal categories, 
based solely on expected health care resource demands. 

For the sixth measure we used a count of the num- 
ber of drugs prescribed to an individual in the period 
between 1 April 2004 and 31 March 2005 as a proxy 
measure of multimorbidity. In an earlier study from 
North America in older adults, a similar measure 
appeared to predict health care utilization better than 
the Charlson index or ACGs. 10 We counted the number 
of unique British National Formulary (BNF) 17 codes 
appearing in the individual's prescription drug data. 
Each code represents one sub-heading within the BNF 
and includes drugs that are in the same class (for exam- 
ple thiazide diuretics have one code and loop diuret- 
ics have another). Therefore, repeated prescriptions of 
the same or very similar medication, including different 
doses or formulations, were only counted once. 

Analyses 

We investigated the explanatory power of each mul- 
timorbidity measure in relation to two outcome 



174 



Family Practice— The International Journal for Research in Primary Care 



measures: (i) number of consultations and (ii) mortality. 
Both outcomes were measured over the 3-year period 
beginning 1 April 2005. Using 3 years of data improved 
the reliability of estimates, particularly for mortality, 
which is a rare event. We used the likelihood ratio test 
(backward selection with a significance level of 0.01) to 
determine which covariates to include in the baseline 
model. We considered age, sex, age-by-sex interactions, 
deprivation and GP practice effects. Age, at 1 April 
2005, was categorized into 10-year age bands with the 
exception of 18-29 and 90+ as lower and upper catego- 
ries, respectively. Deprivation was categorized into 10 
deciles based on the Index of Multiple Deprivation 
(IMD) 2007, which incorporates seven dimensions of 
deprivation and relates to the individual's Lower Layer 
Super Output Area (LSOA). 18 

Each multimorbidity measure was included in the 
model as a categorical variable to allow for possible 
non-linearity in the relationship between multimorbid- 
ity and the relevant outcome. For each measure we col- 
lapsed the uppermost categories to ensure there were 
enough individuals in each cell to allow estimation of 
the parameters. Spearman rank correlations were used 
to assess the level of agreement between the numerical 
measures of multimorbidity used in the models. 

The number of consultations over the 3 years follow- 
ing 1 April 2005 was modelled using Generalised Linear 
Models (GLMs) assuming a log link function and nega- 
tive binomial distribution. We included consultations 
with GPs and practice nurses, both face-to-face and by 
telephone. For any individual who was not observed for 
the entire 3-year period due to death or attrition (e.g. 
transferring out of the practice) we applied a weighting 
to the observed number of consultations (the weighting 
was equal to the inverse of the proportion of the 3-year 
period during which they were observed). 

Mortality was modelled using logistic regression 
models where the outcome was presence, or absence, 
of death during the 3 years following 1 April 2005. 
We considered using a Cox proportional hazards 
regression model; however, the extremely high pro- 
portion of censoring in our data (more than 96% of 
the sample did not die during the three-year period) 
meant that standard measures of model fit, such as 
R-squared based measures, would have been unrelia- 
ble. 19 As age is such a strong predictor of mortality we 
also ran the mortality regression models using a sub- 
set of the sample, namely individuals aged 65 years 
and over. In this model we included age in years as 
a categorical variable (i.e. one coefficient estimated 
for each year rather than each 10-year age band). The 
results from this analysis were then compared with 
the results from the logistic regression model using 
the full sample. 

Models were compared using Akaike's Information 
Criterion (AIC) 20 and Bayesian Information Criterion 
(BIC), 21 for both of which a smaller value provides an 



indication of a better fitting model. The BIC provides a 
greater penalty for additional parameters in the model 
than the AIC. We also calculated a deviance-based 
R-squared measure 22 that may be interpreted as the 
proportion of uncertainty in the relevant outcome that 
has been explained by the model. 

Software 

We used John Hopkins ACG System Version 8.2 to 
obtain the EDC, ACG and RUB classifications 16 and 
STATA Version 11.2 for the statistical analysis. 



Results 

Summary statistics 

Of the 95372 individuals in the original sample, 95 188 
(99.8%) were included in the complete case analysis 
after excluding 184 individuals with missing depriva- 
tion data. The mean age was 49 years (SD 18 years) and 
51% of the sample was female. 

Table 1 shows descriptive statistics for each outcome 
in relation to age, sex and deprivation. The number of 
consultations per patient per annum ranged from 0 to 
113 with mean, SD and median of 4.3, 5.0 and 3, respec- 
tively. Of the 95372 people in the full sample, 3198 
(3.4%) died and 10081 (10.6%) left their practice dur- 
ing the 3-year observation period. 

Figure 1 shows the frequency distribution of each 
measure of multimorbidity (with the exception of the 
ACG categories). The vertical axis shows the percent- 
age of the sample at each level or category of multi- 
morbidity. The QOF disease count and Charlson index 
each show an extremely positively skewed distribution, 
with the majority of people having a count or score of 
zero. The EDC count and number of prescribed drugs 
are also positively skewed, but have a greater number 
of categories and a more even distribution. In terms of 
Resource Utilisation Bands very few people were in 
the high utilisation bands. Of the 68 populated ACG 
categories, the largest 12 categories contained over 
80% of the sample. 

Table 2 provides Spearman rank correlations between 
the numerical measures of multimorbidity used in the 
models. There is stronger agreement between the QOF 
disease count, EDC count and number of prescribed 
drugs (r ranging from 0.60 to 0.65) than there is with 
the Charlson index score (r ranging from 0.42 to 0.49). 

Number of consultations 

Table 3 gives the deviance-based R-squared, AIC and 
BIC values for the different models predicting the 
three-year consultation rate. A model containing age, 
sex, age-by-sex interaction, deprivation and GP prac- 
tice alone explained 22% of the observed uncertainty. 
Inclusion of the number of prescribed drugs, ACG cat- 
egory, EDC count, RUB category, QOF disease count, 
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Table 1 Descriptive statistics on the number of consultations per 
patient per annum, and the number of deaths, over a 3-year period 



Variable Total Number of Number of 

frequency consultations deaths, n (%) 
per patient per 
annum, mean 
(SD) [median] 



Age category 



18-29 


16 099 


2 7 (3 4) \2\ 


27 (0.2%) 


30-39 


17 338 


3 0 (3 7*1 121 


33 (0.2%) 


40-49 


18 059 


3 3 (4 H 121 


89 (0.5%) 


50-59 


16 220 


4 2 (4 5") 131 


232 (1.4%) 


60-69 


12 469 


5 8 (5 5) 151 


411 (3.3%) 


70-79 


9227 


7.5 (6.3) [6] 


823 (8.9%) 


80-89 


4957 


7 6 (7 0) 161 
/.u v'-uy L U J 


1137 (22.9%) 


90+ 


1003 


5 6 (5 91 141 


446 (44.5%) 


Gender 








Male 


46 768 


3.4 (4.6) [2] 


1499 (3.2%) 


Female 


48604 


5.1 (5.2) [4] 


1699 (3.5%) 


Deprivation 








f (least) 


9625 


4.1 (4.6) [3] 


236 (2.5%) 


2 


9417 


4.1 (4.7) [3] 


288 (3.1%) 


3 


9654 


4.2(5.0) [3] 


298 (3.1%) 


4 


9390 


4.1 (4.6) [3] 


298 (3.2%) 


5 


9512 


4.0 (4.8) [3] 


271 (2.8%) 


6 


9537 


4.4 (5.1) [3] 


322 (3.4%) 


7 


9554 


4.4 (5.2) [3] 


341 (3.6%) 


8 


9464 


4.3 (5.0) [3] 


335 (3.5%) 


9 


9541 


4.5 (5.2) [3] 


369 (3.9%) 


10 (most) 


9494 


4.7 (5.5) [3] 


429 (4.5%) 


Overall 


95 372 


4.3 (5.0) [3] 


3198 (3.4%) 



Total number of subjects with non-missing deprivation data is 95 188. 
Data source: GPRD which has now been incorporated into the 
Clinical Practice Research Datalink (http://www.cprd.com). 



or Charlson index score increased this to 42%, 37%, 
36%, 35%, 30% and 26%, respectively. Both AIC and 
BIC showed the same pattern of support for the models 
as the deviance-based R-squared values. 

Mortality 

Table 4 gives the deviance-based R-squared, AIC and 
BIC values for the logistic regression models predict- 
ing the probability of death within 3 years. The model 
containing the ACG categories could not be reliably 
estimated as several of the less populated ACG catego- 
ries had no observed deaths. A model containing age, 
sex and deprivation alone explained 28% of the uncer- 
tainty. The R-squared value suggested there was very 
little improvement from including any of the measures 
of multimorbidity, with gains of no more than 3% over 
the baseline model. Further investigation showed that 
a model with age category alone (R-squared = 27.5%, 
further results not shown) suggested that age was the 
major contributor to any uncertainty explained in the 
outcome. Nonetheless, the deviance-based R-squared, 
AIC and BIC values all showed the Charlson index 



score as the highest ranked model. This was followed 
by a model containing the number of prescribed drugs. 
When considering ranking of the models beyond that 
point, there is no clear distinction in R-squared values, 
nor is there complete agreement between AIC and BIC 
values. 

Fitting the logistic regression models using only indi- 
viduals aged 65 years and over showed similar results. 
The model which included the Charlson index score 
performed best, followed by the number of prescribed 
drugs. Even using this limited age range, improvements 
in the deviance-based R-squared values due to the 
inclusion of a multimorbidity measure were still small 
(ranging between 1.3% and 3.5%). 

Discussion 

Summary of main findings 

This research has demonstrated that all of the multi- 
morbidity measures had moderate predictive validity in 
relation to consultation rates in primary care, with the 
number of prescribed drugs having the greatest predic- 
tive validity followed by the three ACG based meas- 
ures. With regard to mortality, none of the measures 
added greatly to the fit of a model containing demo- 
graphic variables, but the Charlson index performed 
best, followed by the number of prescribed drugs. The 
other diagnosis-based measures of multimorbidity (the 
ACG-related measures and the number of chronic con- 
ditions included in the QOF) added little to the model. 

The Charlson index was originally designed to pre- 
dict mortality in hospitalized patients, and these results 
confirm that it is similarly valid in the community popu- 
lation. It may not be surprising that the number of drugs 
predicts consultation rates, since people who consult 
more often may also be prescribed more drugs, but it is 
notable that this measure also predicts mortality better 
than diagnosis-based measures other than the Charlson 
score. This suggests that this relatively simple measure 
may be a useful proxy measure for primary care-based 
studies that are exploring a range of outcomes and 
need a case-mix measure to account for multimorbidity 
amongst participants. 

Comparison with previous studies 
Our findings build on previous research by enabling the 
comparison of several different approaches to measur- 
ing multimorbidity in the same population, and com- 
paring their predictive validity in relation to both a 
process and an outcome variable. Our research is also 
based on a much larger and more reliable dataset than 
most previous studies. 

Several previous studies have confirmed that the 
Charlson index is a valid measure to predict mortality 
amongst community populations. 91011,23 Some recent 
studies have suggested that the Elixhauser index 24 may 



176 



Family Practice— The International Journal for Research in Primary Care 





o 




CD 


0) 


O 




LO 


S 


o 


sa 


«* 




o 


o 


ro 


0) 






O 


<c 


eg 


c 


o 


d) 




o 




: — 




0) 


o 


Q- 






12 3 4 
QOF chronic disease count 



5+ 





o 


' — ' 








"5. 




— 


o 

CO 


(0 


w 






O 


o 


CM 


0) 




CT> 




OT 


O 






C 




0) 




o 






O - 


0) 




0. 






Non-user Low morbidity High 

Healthy user Moderate Very high 

Resource Utilisation Band 





© 








o 




CD 


£ 


O 




LO 




O 


M- 

o 






o 


0) 
CD 


oo 


TO 


o 




CM 


C 




<D 


O 


O 




i_ 




a> 


O 


a. 






~ I 1 1 1 1 1 — 

0 1 2 3 4 5+ 
Charlson Index score 



a> 
a. 



a> 
o> 
to 

a> 

o 

: 

a. 



o 



O 

C-J 



DDDdd 



t — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



0 2 4 6 
13 5 7 



10 13-14 17-18 21-25 
9 11-12 15-16 19-20 over25 



Number of prescribed drugs 



a; 



a. 



LO 



o 



a> ° 

O) 



c 

o 



□ □ 



2 4 

3 



~ 1 1 1 1 1 1 1 — 

10 12 14 
9 11 13 15+ 



5 7 

EDC count 
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prescribed drugs. Percentage of the sample is shown on the vertical axis. 



be superior to the Charlson index 112526 but this is cal- 
culated using ICD codes and algorithms have not yet 



Table 2 Spearman rank correlations for the numerical measures of 
multimorbidity 





QOF disease 


Charlson index 


Count of EDCs 




count 


score 




QOF disease count 


1.00 






Charlson index 


0.46 


1.00 




score 








Count of EDCs 


0.65 


0.49 


1.00 


Number of 


0.60 


0.42 


0.65 


prescribed drugs 









been developed for use with Read codes, so it was not 
included in this study. 

Our findings are compatible with previous research 
in a range of populations on prediction of different 
aspects of health care utilization or health care costs. 
Both simple counts of number of distinct medications 
in the previous year and more sophisticated medica- 
tion-based measures such as the Chronic Disease Score 
or RxRisk outperform the Charlson index in predicting 
health care utilization. 1011,27 ' 28 . 

Strengths and limitations 

The GPRD provides a large, reliable and validated data- 
set of patient records in primary care, for a population 
that is broadly representative of England. Most previous 
research has been conducted in administrative databases 



Comparing measures of multimorbidity in primary care 



177 



Table 3 Model fit statistics for negative binomial regression models 
predicting 3-year consultation rate 



Model 



Deviance-based R 2 



AIC 



BIC 



No multimorbidity 21.9% (7) 
measure 

QOF disease count 30.2% (5) 

Charlson index 25.9% (6) 
score 

Count of EDCs 36.1% (3) 

ACG category 37.3% (2) 

Resource 34.7% (4) 
Utilisation Band 

Number of prescribed 41.6% (1) 
drugs 



648367 (7) 650373 (7) 



636666 (5) 
642855 (6) 

627617 (3) 
625798 (2) 
629809 (4) 



638720 (5) 
644908 (6) 

629766 (3) 
628438 (2) 
631863 (4) 



618632(1) 620799 (1) 



All models include age category, sex, age-by-sex interaction, depriva- 
tion (based on IMD2007) and GP practice as fixed effects. 
The associated ranking of each model is shown in parentheses, with 1 
representing the best performing model and 7 representing the worst 
performing model. 

For AIC and BIC a smaller value indicates better model fit. 



Table 4 Model fit statistics for logistic regression models predicting 
the probability of death over a 3-year period 



Model 



Deviance- 
based R 2 



AIC 



BIC 



No multimorbidity 


28.2% 


(6) 


20079 


(6) 


20250 


(6) 


measure 














QOF disease count 


29.8% 


(3) 


19636 


(3) 


19854 


(3) 


Charlson index score 


31.3% 


(1) 


19226 


(1) 


19443 


(1) 


Count of EDCs 


29.8% 


(3) 


19667 


(4) 


19979 


(5) 


Resource Utilisation Band 


29.5% 


(5) 


19728 


(5) 


19946 


(4) 


Number of prescribed drugs 


30.9% 


(2) 


19361 


(2) 


19693 


(2) 



All models include age category, sex, and deprivation (based on 
IMD2007) as fixed effects. 

The associated ranking of each model is shown in parentheses, with 1 
representing the best performing model and 7 representing the worst 
performing model. 

For AIC and BIC a smaller value indicates better model fit. 



of uncertain reliability and/or on specific sub-popula- 
tions. The findings from this study are therefore likely to 
provide a robust indication of the validity of these meas- 
ures and be applicable to other primary care research in 
England and possibly in other developed countries. 

There may be limitations due to how the measures 
were constructed. Both the Charlson index and the 
ACG measures are based on ICD codes, and although 
algorithms have been developed to calculate them from 
Read coded data, few previous studies have used them 
in this way. The findings may also be sensitive to deci- 
sions made in operationalizing measures. For example, 
we calculated a medication measure based on the num- 
ber of prescriptions for drugs with different BNF codes, 
while other researchers have counted the number of 



prescription items. 10 Although the number of different 
drugs (rather than the number of prescription items) 
may be a better proxy for the number of conditions 
an individual has been diagnosed with, the number of 
repeated prescriptions of a drug may be more likely 
to relate to chronic conditions. Similarly, the ACG 
software can be used in several ways. For example, we 
could have included diagnoses recorded over the pre- 
vious 3 years rather than 1 year in calculating ACG 
categories or could have based our analyses on ADGs 
rather than ACGs. Therefore, it may be worthwhile in 
future studies to investigate the relative performance of 
different approaches to counting prescriptions or using 
the ACG software. Finally it should be noted that all 
measures of multimorbidity (and all epidemiological 
studies) based on routine data are reliant on whether 
diseases have been diagnosed and coded. 

It is notable that none of the measures used in this 
study added very substantially to age alone in predict- 
ing mortality over the subsequent 3 years. Our find- 
ings are likely to be a consequence of the fact that in 
the population as a whole mortality is a rare event, 
and this may partly explain the contrast between our 
findings and previous research, which has been based 
on elderly populations. 91023 However, we found simi- 
lar results using a subset of the sample aged 65 years 
and over. Although multimorbidity measures are com- 
monly used to adjust for case-mix in epidemiological 
and observational studies, these findings highlight 
that they have only a limited ability to control for 
confounding. 27 

It is notable that the number of prescribed drugs 
over the previous year is not only the most powerful 
predictor of primary care consultations over the sub- 
sequent 3 years but is also (after the Charlson index) 
the second most powerful predictor of mortality. This 
is particularly interesting given that the Charlson index 
was specifically designed to predict mortality, yet has 
limited additional predictive validity over a simple 
count of drugs prescribed. The latter measure is simple 
to use and does not require statistical expertise or pro- 
prietary software to calculate. Future research should 
explore whether this measure predicts other outcomes 
such as quality of life, self-rated health and functional 
ability. 



Conclusion 

No single measure adds considerably to age alone in 
predicting mortality amongst the adult population 
consulting in primary care. All the measures included 
in this study have weak predictive validity for this pur- 
pose but the Charlson index is the most powerful. The 
number of prescribed drugs is the most powerful meas- 
ure for predicting future health care consultations and 
is the second most powerful measure for predicting 
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mortality. It may have potential as a useful and simple 
proxy measure of multimorbidity in primary care. 
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